Magnetic fields generated by r-modes in accreting 

O ■ millisecond pulsars 
(N 

■ Carmine Cuofano and Alessandro Drago 

Dipartimento di Fisica, Universita di Ferrara 
and INFN sez. Ferrara, 44100 Ferrara, Italy 

CO ' E-mail: cuofano@fe.infn.it 

Abstract. In millisecond pulsars the existence of the Coriolis force allows the development 

■ of the so-called Rossby oscillations (r-modes) which are know to be unstable to emission of 
pl"^ ' gravitational waves. These instabilities are mainly damped by the viscosity of the star or by the 

• , existence of a strong magnetic field. A fraction of the observed millisecond pulsars are known to 

i-^ ■ be inside Low Mass X-ray Binaries (LMXBs) , systems in which a neutron star (or a black hole) 

is accreting from a donor whose mass is smaller than 1 Mq. Here we show that the r-mode 
^ , instabilities can generate strong toroidal magnetic fields by inducing differential rotation. In 

i-^ ' this way we also provide an alternative scenario for the origin of the magnetars. 

^ ' 1. Introduction 

0^ , The r-mode oscillations in all rotating stars are unstable for emission of gravitational waves |T] . 

These modes play therefore a very important role in the astrophysics of compact stars and in the 

■ search for gravitational waves. On the other hand the existence of millisecond pulsars implies 
the presence of damping mechanisms of the r-modes. Damping mechanisms are associated with 

' bulk and shear viscosity and with the possible existence of the so called Ekman layer. The latter 

0^ . is located at the interface between the solid crust and the fluid of the inner core and in this 

region friction is significantely enhanced respect to friction in a purely fluid component. All 
these mechanisms are strongly temperature dependent. 
^ , An important class of rapidly rotating neutron stars are the accreting millisecond pulsars 

\—i I associated with Low Mass X-ray Binaries (LMXBs). For these objects the internal temperature 

is estimated to be in the range IqS-IO^-^ K [a |3] and their frequencies can be as large as ~ 650 
Hz. In this range of temperatures and in the case of a purely nucleonic star, bulk and shear 
viscosities alone cannot stabilize stars whose frequency exceeds ~ 100 Hz. A possible explanation 
of the stability of stars rotating at higher frequencies is based on the Ekman layer, but recent 
calculations show that this explanation holds only for rather extreme values of the parameters 
[4j. In this contribution we propose a new damping mechanism based on the generation inside 
the star of strong magnetic fields produced by r-mode instabilities. This same mechanism has 
been proposed in the case of rapidly rotating, isolated and newly born neutron stars in jFj "B]. 
In that paper the mechanism which generates the magnetic field is investigated only during the 
relatively short period in which the star remains always in the instability region. In our work 
we consider accreting stars and we investigate the interplay between r-modes and magnetic field 
on an extremely long period and we show that in this scenario a very strong magnetic field can 



be produced. 



2. R-mode equations in the presence of magnetic field 

R-mode instabilities are associated to kinematical secular effects wich generate differential 
rotation in the star and large scale mass drifts, particularly in the azimuthal direction. 
Differential rotation in turn can produce very strong toroidal magnetic fields in the nucleus 
and these fields damp the instabilities extracting angular momentum from the modes. In order 
to derive the equations regulating the evolution of r-modes in the presence of a pre-existent 
poloidal magnetic field we have modified the equations derived in |7j, taking into account also 
the magnetic damping. 

We use the estimate given in p] for the gravitational radiation reaction rate due to the I = m = 2 
current multipole 

Fg = ^M,,,Rf,Pll s (1) 
as well as the bulk and shear viscosity damping rates 

s 

(2) 

where we have used the notation -/Vfi.4 = M/IAMq, Riq = R/10 Km, P^-^ = P/1 ms and 
Tg = T/IO^ K. 

The total angular momentum J of a star can be decomposed into a equilibrium angular 
momentum J* and a perturbation proportional to the canonical angular momentum of a r-mode 
Jc-. 

J = J^{M, + Kj)Jc, Je = -K^a^,h (3) 
where K(^j f.-^ are dimensionless constants and J* = I^Vt. 

Following Ref.[8] the canonical angular momentum obeys the following equation: 

dJjdt = 2J,{Fg{M, n) - [F,{M, n, T,) + Fm^ (M, Q, Bp)]} (4) 

where F^ = Fg + Fb is the viscous damping rate and we have introduced the magnetic damping 

rate Frm that we discuss in the next section. 

The total angular momentum satisfies instead the equation: 

dJ/dt = 2J,Fg + ja{t) - hnPm^ (5) 

where Ja is the rate of accretion of angular momentum and we have assumed it to be 
Ja = M{GMRy^'^ , and Fm^ is the magnetic braking rate associated to the poloidal magnetic 
field. 

Combining the equations ^ and ([5j) than we give the dynamical evolution relation 

^ = a{Fg-F^-F^J+a[K,Fg + {l-K,){F^ + F^J]K,a'-^(^j^y\^ 

^ = -2K,na'[K,Fg + {l-K,m + F^J]-^ + ^i^-^Y' -nF^^ (6) 

where I = IMR^ with I = 0.261 for an n=l polytrope. The previous equations can be simplified 
when the star is close to the instability region identified by the condition Fg — F^ — Fm ~ and 
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we obtain: 



dt ''^ ™^ M / \MR^J 

dot „ ,^ ^ ^ , 2 N olM ( G \ 

We note that ii'c = 9.4 x 10~^ and that in this approximation the value of Kj is unimportant. 
Moreover the magnetic breaking is negUgible for magnetic fields typical of the LMBXs (lO'^- 
lO^G) and for accretion rates considered (10~'^-10~^ M0yr~^). 

Viscosity depends critically on temperature. We include three factors in modelling the 
temperature evolution: modified URCA cooling, shear viscosity reheating, and accretion heating. 
The cooling rate due to the modified URCA reactions, e^, is given in [9] 

= 7.5 X 10=^9 Mi/4^r| erg s'^ (8) 

The neutron star will be heated by the action of shear viscosity on the r-mode oscillations. The 
heating rate due to shear viscosity, e^, is given by [Ij 

e, = 2a^n^MR^JFs = 8.3 x lO^'^a'^Q'^JM^f^R'^^^^Tg'^ erg s"^ (9) 

where J = 1.635 x 10"^ 

Accretion heating have two components. We use the estimates given in [10]. The first 
contribution arises when accreting matter undergoes nuclear burning at the surface of the star 

e„ = — X 1.5MeV = 4 x lO^^Mi 4 erg s"^ (10) 

ruB 

where is the mass of a barion. 

The second contribution arise becouse the flow is assumed to be advection dominated. Matter 
falling in towards a star liberate ~ GM/R of potential energy per unit mass. In a non-advection 
dominated flow most of this would be dissipated as heat in accretion disk, while in an advection 
dominated flow this energy is carried in with the flow of matter. The heating rate is 

RGMM .1 13/6 -1 

eh ~ ^ = 8 X 105^MiTMi.4 ergs-i (11) 

Finally we use the estimate of the heat capacity given in [TU] 

a = 1.6 X lO^^Ml^^Tg erg K"^ (12) 



3. Magnetic damping 

The crucial ingredient introduced in the previous Section is the magnetic damping rate, which we 
have inserted in the evolution of the r-modes. This specific modification of the r-mode equations 
is at the basis of the phenomenolgy we are going to discuss and was never considered in previous 
calculations. The expression of the magnetic damping rate has been derived in [5, i6j, where it 
has been shown that while the star remains in the instability region, the r-modes generate a 
differential rotation which can greatly amplify a pre-existent magnetic field. More specifically, 
if a poloidal magnetic field was originaly present, a strong toroidal field is generated inside the 
star. The energy of the modes is therefore transfered to the magnetic field and the instability is 
damped. 



The expression of the magnetic damping rate reads: 




4(1 -p) 



(13) 



9ttp ■ (8.2 X 10-3) 



where p and A' are dimensionless parameters of order unit. The time integral over the r-mode 
amphtude a takes contribution from the period during which the star is inside the instabihty 
region. The crucial point is if the toroidal magnetic field generated during an instability phase 
can be stabilized or unwind on a Alfven timescale. 



We consider a scenario in which the mass accretion spin up a initially slowly rotating neutron 
star and we investigate the case in which the star cannot find a state of thermal equilibrium, 
then it enters a regime of thermogravitational runaway within a few hundred years of crossing 
the r-mode stability curve. In Figure ^ and Figure ([2]) we show the evolution of temperature, 
spin frequency and r-mode amplitude obtained without magnetic field. At the instant A the 
accreting neutron star enters the r-modes instability window. Until the star is close to equilib- 
rium the instability grows slowly and the star moves to higher temperatures mainly because of 
non thermal equilibrium. After the instant B the r-mode instability grows exponentially due to 
the decrease of the shear viscosity with increasing temperature. As a consequence the r-modes 
amplitude rapidly reaches the saturation value and the viscosity heats significantly the star. At 
this stage the star loses angular momentum by emission of gravitational waves and within few 
hundred years goes out of the instability region (instant C). 

Taking into account the magnetic field, the evolutionary scenario for the star is quite different. 
After the instant B, the secular effects are extremely large and the toroidal magnetic field is 
either produced or amplified by the wrapping of the poloidal magnetic field produced by the 
(mostly) toroidal secular velocity field. Our results obtained solving Eqs.d?]) and reported in 
Figure ([3]) show how the new magnetic field so generated reaches a value Bfor ~ 10^^ G in about 
one hundred years and stabilizes the star when it is still in the instability region. 
The subsequent evolution of the system is difficult to control because it strongly depends on the 
evolution of the magnetic field. Analytical and numerical simulations show that strong purely 
toroidal fields are unstable. We will limit ourselves to outline a possible scenario that we dis- 
cuss in a forthcoming paper. Following the results of some numerical simulations, the magnetic 
field can evolve on an Alfven timescale into a stable configuration of a mixed poloidal-toroidal 
twisted-torus shape |llj . In this scenario the newly formed magnetic field can provide a effective 
damping of the r-modes and rapidly rotating stars can be stabilized also in the absence of the 
Ekman layer. 



4. Results 



It is a pleasure to thank Giuseppe Pagliara and Andreas Reisenegger for many useful and 
stimulating discussions. 




Figure 1. Path followed by the accreting neutron star, without toroidal magnetic field, in 
the Temperature vs Frequency plane. Here M = 10~^ Mq yr~^ and poloidal magnetic field 
Bpoi = G. 
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Figure 2. Temporal evolution, without toroidal magnetic field, of frequency, r-mode amplitude 
and temperature of the accreting neutron star. 
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Figure 3. Temporal evolution, with toroidal magnetic field, of frequency, r-mode amplitude 
and temperature of the accreting neutron star. 
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